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Overview 


This report covers the research work performed under Grant 
NAG- 1-637 for the period ending February 28, 1989. A detailed de- 
scription of the fracture analysis of transverse crack tip delaminations 
is presented in the following sections. This work was performed 
during the first six months of the grant period and has been accepted 
for presentation at the 30th Structures, Structural Dynamics and 
Materials Conference (Mobile, Alabama, April 1989). The following 
sections are adapted from the aforementioned paper. 

Abstract 

Delamination is a predominant failure mode in continuous fiber 
reinforced laminated composite structures. One type of delamination 
is the transverse crack tip delamination which originates at the tip of 
transverse matrix cracks. An analytical model based on the sublami- 
nate approach and fracture mechanics is developed in this paper to 
study the growth of such delaminations. Plane strain conditions are 
assumed and estimates are provided for the total strain energy release 
rate as well as the mode I and mode II contributions. The energy re- 
lease rate estimates are used in combination with a simple failure law 
to predict critical delamination growth strains and stresses. These 
predictions are compared with experimental data on T300/934 
Graphite Epoxy [±25/90 n ] s laminates in the range n=.5 to 8. A good 
agreement is demonstrated for the range of n where the experimental 
observations indicate transverse crack tip delamination to be the pre- 
dominant failure mode. 
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Introduction 


Fiber reinforced composites are now being used in a wide variety 
of engineering structures. The concept of directional strength and 
stiffness has been, for the most part, understood sufficiently to enable 
efficient load bearing designs. One of the current major issues in 
composite structures is the understanding and prediction of damage 
modes and failure mechanisms. A thorough knowledge of the failure 
mechanisms is bound to lead to the design of efficient and durable 
structures. Failures in these materials often initiate in the form of ma- 
trix cracks or delaminations. Matrix cracks refer to intralaminar fail- 
ures whereas delaminations refer to interlaminar failures. 

Matrix cracks usually occur within laminates where the fibers 
run at an angle to the primary load direction. Hence, such matrix 
cracks are also called transverse cracks. Based on the location and di- 
rection of growth, two distinct types of delamination can be discerned. 
These two types are called edge delamination and local or transverse 
crack tip delamination. Edge delaminations initiate at the load free 
edges of the laminate whereas local delaminations start from a trans- 
verse matrix crack. In many cases, both types occur concurrently with 
varying levels of interaction. It has been observed in simple tension 
tests of uniform rectangular cross section specimen (Edge 
Delamination tests) that delaminations initiate along the load free 
edges and propagate normal to the load direction. Transverse matrix 
cracks running parallel to the fibers have also been observed in off axis 
plies such as 90° plies. Such transverse cracks terminate where the 
ply orientation changes. Delaminations can also originate at the inter- 
faces where transverse cracks terminate. These delaminations, called 
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transverse crack tip delaminations or local delaminations, grow nor- 
mal to the transverse crack from which they originate. In the case of 
90° plies, the growth direction is parallel to the load. 

The growth process of edge delaminations and local delamina- 
tions is often modelled using a fracture mechanics approach leading to 
the calculation of a strain energy release rate. This is because the 
strain energy release rate can correlate delamination behavior from 
different loading conditions and can account for geometric depen- 
dencies. The strain energy release rate associated with a particular 
growth configuration is a measure of the driving force behind that fail- 
ure mode. In combination with appropriate failure criteria, the strain 
energy release rate provides a means of predicting the failure loads of 
the structure. 

Several methods are available in the literature for analyzing edge 
delaminations. These include finite element modelling 1 * 3 , complex 
variable stress potential approach 4 , simple classical laminate theory 
based technique 5 and higher order laminate theory including shear 
deformations 6 . Finite element models provide accurate solutions but 
involve intensive computational effort. Classical laminate theory (CLT) 
provides simple closed form solutions and is thus well suited for pre- 
liminary design evaluation. However, CLT provides only the total en- 
ergy release rate, and thus, in a mixed mode situation, there is insuf- 
ficient information to completely assess the delamination growth ten- 
dency. A higher order laminate theory including shear deformations 
has the ability to provide the individual contributions of the three frac- 
ture modes while retaining the simplicity of a closed form solution. A 
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shear deformation model is available for edge delamination and has 
been shown to agree well with finite element predictions 6 . 

Crossman and Wang 7 have tested T300/934 Graphite epoxy 
[±25/90 n ls specimens in simple tension and reported a range of be- 
havior including transverse cracking, edge delamination and local de- 
lamination. O'Brien 8 has presented classical laminate theory solutions 
for these specimen, demonstrating reasonable agreement in the case 
of edge delamination but with some discrepancies in the local delami- 
nation predictions. A finite element model combining edge and local 
delaminations has been proposed by Law 9 . His predictions, however, 
do not fully explain the dependency of the critical strain on the num- 
ber of 90° plies. 

In this paper, a shear deformation model is developed for the 
analysis of local delaminations originating from transverse cracks in 
90° plies located in and around the specimen midplane. Plane strain 
conditions are assumed and thickness strain is neglected. Delamina- 
tions are assumed to grow from both ends of the transverse crack tip. 
The transverse crack is treated as a free boundary and the delamina- 
tion is considered to be the crack whose growth behavior is to be 
modelled. The sublaminate approach 10 - 11 is used to model different 
regions of the specimen. The resulting boundary value problem is 
solved to obtain the interlaminar stresses, total strain energy release 
rate and energy release rate components. Critical local delamination 
growth loads are predicted for [±25/90 n ls specimens. 

Analytical Model 

The formulation is based on the sublaminate approach detailed 
in Ref. 10. A longitudinal section illustrating the geometry of a generic 
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configuration is shown in fig. 1. The central region is assumed to be 
made of 90° plies with an isolated transverse crack in the middle. 
Delaminations are assumed to grow from both ends of the transverse 
crack, and towards both ends as shown. From symmetry considera- 
tions, only one quarter of the configuration is modelled. The modelled 
portion is divided into four sublaminates as shown in fig. 2. The top 
surface (sublaminates 1 and 4) is stress free. In order to simplify the 
analysis, plane strain conditions are assumed and the thickness strain 
(£ z ) is set to zero. The consequence of this, combined with the fact 
that the w displacement is zero along the center line, is that w is zero 
in sublaminates 1, 2 and 3. Also, this approximation does not allow for 
the enforcement of boundary conditions on the shear stress resultants, 
leading to incorrect estimates of the interlaminar normal stresses. 
The interlaminar shear stresses, however, are not affected by this as- 
sumption 6 * 10 . These assumptions lead to considerable simplifications 
in the analysis. In spite of the simplifications, reliable energy release 
rate components can be estimated based on the interlaminar shear 
stress distributions 6 * 10 . 

A generic sublaminate is shown in fig. 3 along with the notations 
and sign conventions. The peel and interlaminar shear stresses are 
denoted by P and T, respectively, with t and b subscripts for the top 
and bottom surfaces, respectively. The axial stress resultant, shear 
stress resultant and bending moment resultant are denoted by N, Q 
and M, respectively. A summary of the governing equations is pre- 
sented in the following paragraphs for convenience. These equations 
are derived for a generic sublaminate using the principle of virtual 
work in Ref. 12. 
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The x and z displacements within the sublaminate are assumed 
to be of the form 

u(x,z) = U(x)+zp(z) (1) 

w(x,z) = W(x) (2) 

Here U represents the axial midplane stretching and W is the trans- 
verse displacement. The shear deformation is recognized through the 
rotation (3. The origin of the coordinate axes for the sublaminates is 
taken at the delamination tip as shown in fig. 4. The equilibrium 
equations take the form 

N, x +Tt-Tb =0 (3) 

Qoc+PfPb = 0 (4) 

M, x -9+|(Tt+Tb) = 0 (5) 

where h is the thickness of the sublaminate. The constitutive rela- 
tions in terms of the force and moment resultants are 

N = AiiU,x+Bii|3,x (6) 

Q=A5 5 (P+W,x) (7) 

M = BiiU.x+DnP.x (8) 

where Ay, By and Dy are the classical laminate theory axial, coupling 
and bending stiffnesses, respectively. The boundary variables to be 
prescribed at the sublaminate edges are 

N or U 
M or P 
Q orW 

Additionally, at the interfaces between sublaminates, reciprocal 
traction and displacement matching boundary conditions have to 
specified. 
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Solution Procedure 


A detailed solution is provided in the Appendix. A summary is 
provided in this section for convenience. The variables are subscripted 
to indicate the sublaminate in which they occur. The solutions in 
sublaminates 1 and 2 are coupled by the reciprocal interlaminar 
stresses denoted Ti and Pi and by displacement continuity at the 
common interface. Assuming exponential solutions for the axial force 
and bending moment resultants leads to an eigenvalue problem involv- 
ing the exponential parameter s. The eigenvalues turn out to be 0 and 
two nonzero values (say si and S2) occurring in positive and negative 
pairs. Since the response decays from the delamination (crack) tip, 
only the exponentially decaying terms are considered in the solutions. 

The following boundary conditions from the ends of the mod- 
elled region are enforced. 


N 2 ( 0 ) = 0 ( 9 ) 

94(a) =0 (10) 

Ma)=0 (11) 

N 1 +N 2 = Applied Load (12) 


Further, the following displacement matching conditions are 
applied. 


(,^) = u 2 (,£) 

(13) 

Ui(0) = U 4 (0) 

(14) 

u 2 ( 0 ) = u 3 ( 0 ) 

(15) 

(3 1(0) = WO) 

(16) 


It should be noted that a $2 and P3 matching condition cannot be 
applied at this level of modeling since it would amount to specifying 
both W and Q 6 - 12 . Consequently, there is a displacement discontinuity 
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at the delamination tip. The effect of this will be discussed subse- 
quently. To eliminate rigid body displacements, Ui is set to zero at 
the left end. The following solutions can then be obtained for the re- 


sultants in sublaminates 1 and 2 . 

Ni = aie s i x +a2e s 2 x +£Aii(i) (17) 

N 2 = -aie s i x -a2e s 2 x +£Ai i(2) ( 1 8) 

Mi = aikie s i x +a2k2e s 2 x ( 1 9) 

M 2 = aik3e s i x +a2k4e s 2 x (20) 

The interlaminar shear and peel stresses between sublaminates 1 and 
2 can be obtained from equilibrium as 

Ti=aiSie s i x +a2S2e s 2 x (21) 

Pi^ki+^^ajs^esix) 

+ ( k2+ ¥)( a 2 s 2 2eS2X ) (22) 


In the above solutions, the ki parameters are dependent on the 
eigenvalues and the stiffness of sublaminates 1 and 2 , the aj parame- 
ters depend on the ki parameters and the crack length a, and £ is de- 
fined as 

P 1 

e = 2 b Au(i)+Aii( 2 ) (23) 

where P is the uniform axial force applied on the specimen and b is 
the specimen width. Expressions for the eigenvalues and the aj and kj 
parameters are provided in the Appendix. 

Proceeding on to sublaminates 3 and 4, the following solutions 
can be written. 


$ 

It 

0 

(24) 

M 3 = <pi sinh(G)3x)+92 cosh(oc> 3 x) 

(25) 

<P 2 = aik3+a2k4. 

(26) 


where 
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and 


(pi = -q>2 coth(C03a), 

■4 


©3 


f A55(2) 

Dh(2) 


( 27 ) 

( 28 ) 

( 29 ) 

( 30 ) 


N4 = E (Aii(i)+Ah(2)) 

M4 = aiki+a2k2 

The corresponding displacement solutions are provided in the 
Appendix. The compliance of the specimen can be evaluated as 

C = 2^ ,31) 

The total energy release rate, Gt. per crack is then given by 

r H dC f o 91 

° T= 2 b da (32) 

Use of the previously described solutions leads to the following ex- 
pression. 

p 2 


Gt 


_ p2 ( 1 
” 2 P(Ah| 


+I1-I2) 


( 33 ) 


lAu(l) An(i)+Aii(2) 
where the quantities Ii and I2 contain exponential terms dependent 
on the delamination length. Using the virtual crack closure technique, 
from the relative displacements in the cracked portion and the inter- 
laminar stresses ahead of the crack tip, the mode I and mode II en- 
ergy release rate contributions can be obtained. The mode III energy 
release rate is zero from the assumption of plane strain. The mode II 
energy release rate is given by 

, 8 
lim 1 


lira 1 r 

Gil = g _^ 0 — |Ti(x- 8 )Au(x)dx 


( 34 ) 


where 8 is the virtual crack step size. The result of the limiting pro- 
cess is zero if there is no singularity in the stress field 10 . So, the limit 
is usually taken as the crack step size d tends to a small value, say A, 
based on the decay length or the length required to capture the 
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essential features of the stress and displacement fields near the crack 
tip. The decay length is dependent on the eigen values si and S 2 . In 
this study, the value of A has been set to 



since it reasonably fulfills the criterion given above. In a similar fash- 
ion, the mode I energy release rate can be obtained based on the nor- 
mal stress (P) and the w displacements near the crack front. The 
normal (peel) stress estimate is inaccurate due to the absence of 
thickness strain. Hence, an alternate approach was used to estimate 
Gi, the mode I energy release rate. The total energy release rate for 
this problem is made up entirely of Gi and Gn (Gm=0). From an esti- 
mate of Gt and Gn, an estimate for Gi can be obtained simply as 

Gi = Gt-Gii (36) 

The critical load for a given specimen can then be evaluated 
based on an appropriate fracture law. This is illustrated in the follow- 
ing section. 

Results and Discussion 

The solutions derived in the previous section have been used to 
model the behavior of [±25/90 n ]s T300/934 Graphite Epoxy specimen 
for n values of .5,l,2,3,4,6,and 8. These correspond to the specimen 
tested by Crossman and Wang 7 . The specimen width and length were 
fixed at .0381m and .015m, respectively, as in the tests and the ap- 
plied uniform axial stress was lOOMPa. The solutions were generated 
using a simple computer program based on the closed form expres- 
sions for the interlaminar stress and energy release rates. 
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An example of the total energy release rate variation with the 
crack length is presented in fig. 5. The asymptotic value of Gt is de- 
noted by Gro hi the figure. It can be observed that after a certain 
crack length, the Gt is independent of the crack length. On the basis 
of curves like the one shown in fig. 5, the crack length was fixed at 10 
ply thicknesses for the remainder of the study. The dependence of 
the mode II contribution of the energy release rate on crack length (a) 
is depicted in fig. 6. Typical interlaminar shear and normal stress 
profiles are presented in figs. 7 and 8, respectively. The correspond- 
ing energy release rates have also been calculated and are presented in 
Table I and fig. 9. 

In order to evaluate the critical loads, an appropriate mixed 
mode fracture law has to be applied, based on the calculated energy 
release components. Since the calculated mode split shows only a 
small variation with n, the simple Griffith criterion Gt=Gtc has been 
used to scale the stresses to obtain the critical delamination growth 
stress (C c ) and strain (e c ) values. The critical energy release rate Gtc 
was chosen as 415 J/m 2 to obtain the critical stresses and strains 
listed in Table I. This value of Gtc is larger than Gi c to account for the 
presence of mode II and the fact that for the material system under 
consideration, Gn c is about four times Gi c . The critical strains are 
plotted against n, the number of 90° plies in fig. 10. The experimental 
results of Ref. 7 and the predictions of Refs. 8 and 9 are also pre- 
sented in the figure for comparison. The predictions of the model de- 
veloped in this paper are represented by the solid line while the ex- 
perimental results are shown as filled squares. The classical laminate 
theory and finite element critical strain predictions of Refs. 8 and 9 
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are represented by triangles with a connecting line and a dotted line 
respectively. The CLT based model agrees well with the shear defor- 
mation model in terms of the total energy release rate. However, the 
CLT based model does not provide information on the mode split and 
thus, the value of G c (=Gi c ) used leads to bias in the predictions. 

In the experiments, the local delamination phenomenon was ob- 
served as the predominant failure mode only for the n=4, 6 and 8 
specimens. The shear deformation model presented in this paper 
provides good agreement with the experimental data in this range. 
For n<4, edge delamination either in the mid plane or in the 25/90 
interface was observed in the tests. Hence, the predictions of the lo- 
cal delamination models in this region are not of consequence as long 
as they do not predict critical loads lower than those predicted by 
edge delamination models. Thus, it can be seen that the shear defor- 
mation model predicts the observed behavior with reasonable accuracy 
and can be used in conjunction with an appropriate edge delamination 
model to predict critical loads accurately for the complete range of n 
values. The edge delamination model presented in Refs. 6 and 12 can 
be used for this purpose. However, a separate model is required to 
account for the midplane (Mode I) edge delamination behavior. The 
development of such a model is described in Ref. 13. 

Conclusions 

A shear deformation model has been developed to analyze local 
delaminations growing from transverse cracks in 90° plies located 
around the mid plane of symmetric laminates. The total energy re- 
lease rate calculations yield the same results as in the case of CLT 
based models. The predictions of the shear deformation model agree 
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reasonably with critical strain experimental data from [±25/90 n ls 
T300/934 Graphite Epoxy laminates. The predicted behavior is such 
that, in combination with an edge delamination model, the critical 
loads can be predicted accurately in the range of n from .5 to 8. 
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Appendix 

Snhlflmlnflte Analysis for Local Dfttaminfltlniifl 

A generic sublaminate is shown in fig. 3 along with the notations 
and sign conventions. The interlaminar normal (peel) and shear 
stresses are denoted by P and T respectively with the t and b sub- 
scripts for the top and bottom surfaces respectively. The axial force 
resultant, shear force resultant and bending moment resultant are de- 
noted by N, Q and M respectively. Plane strain conditions are assumed 
to prevail in the x-z plane and the thickness strain e zz is neglected. 
These assumptions lead to considerable simplification in the analysis. 
The displacements in the x and z directions are assumed to be of the 
form 

u = U(x)+zb(x) (A. 1) 

w = W(x) (A. 2) 

Here U represents the axial stretching and W is the transverse 
(thickness direction) displacement. This formulation recognizes 
shear deformation through the rotation p. The equilibrium equations 
take the form 

N,x+Tt-Tb = 0 (A. 3) 
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Q.x+PfPb = 0 (A. 4) 

M.x-Q+I flt+Tb) = 0 (A. 5) 

where h is the thickness of the sublaminate. The constitutive equa- 
tions in terms of the force and moment resultants are 

N = AnU.x+Bnp,x (A. 6) 

Q=A 55 (P+W,x) (A. 7 ) 

M = BnU,x+Dnp,x (A. 8) 


where A,B and D are the classical laminate theory axial, coupling and 

bending stiffnesses defined in the customary manner as 

h/2 

(Ay.By.Dy) = jCy(l,z,z 2 )dz 
-h/2 

Here, the Cys are the material moduli. For the case of plane strain in 
the x-z plane, the Cs are defined as follows. 

^xx ’Cn C13 0 ' £xx 

• Qzzf = C13C22 0 Itzzf (A. 9 ) 

.TxzJ 0 0 C 55-l Iyxz. 

The boundary quantities to be prescribed at the sublaminate 
edges are 

N or U 
M or p 
Q or W 

Further, at the interfaces between sublaminates, reciprocity of trac- 
tions and continuity of displacements have to be enforced. 

The four sublaminates along with the loads acting on each are 
shown in fig. 4. Setting Pi and Ti as shown automatically satisfies the 
traction matching boundary condition at the 1-2 interface. From 
symmetry, we get w=0 and zero shear stress along the bottom faces of 
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sublaminates 2 and 3. This leads to w=0 in sublaminates 1, 2 and 3. 
Thus, W has been prescribed in these sublaminates and the vertical 
shear force resultant Q cannot be prescribed on these sublaminates. 
Consequently, the calculated peel stress distribution will not be cor- 
rect. In addition, at the 2-3 interface, the (5s cannot be matched, 
since in these sublaminates, specifying p is equivalent to specifying Q 
(through Eq. A. 7). Inspite of these simplifications, reliable energy re- 
lease rate components can be estimated based on the interlaminar 
shear stress distributions. The mode I contribution can then be evalu- 
ated using the total energy release rate, which is not affected signifi- 
cantly by these simplifications. 

For the (±25/90 n )s laminates under consideration. Bn is zero in 
all the four sublaminates. For sublaminates 1 and 2, the equilibrium 
equations and constitutive relationships can be written as 


Ni, x -Ti = 0 

(A. 10) 

N 2 .X+T 1 = 0 

(A. 11) 

Ql,x-Pl = 0 

(A. 12) 

Q2,x+Pl'P2 - 0 

(A. 13) 

Mi )X +TjrTi-Qi = 0 

(A. 14) 

M 2 ,x+^Ti-9 2 = 0 

(A. 15) 

Ni = Aii(i)Ui, x 

(A. 16) 

N 2 = An(2)U2,x 

(A. 17) 

Ql = A55(i)Pi 

(A. 18) 

92 = A55(2)P2 

(A. 19) 

Mi = Dn(i)pi iX 

(A. 20) 

M 2 = Dn(2)p2,x 

(A.21) 
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The subscripts in parentheses refer to the sublaminates to which the 
stiffness coefficients correspond. Eqs. A. 14, A. 15 and A. 12 can be 
rewritten in a modified form as 


hi 

+ ~2 Nl - X = A 55(l)Pl 

(A. 22) 

h2 

M2,x - ~2 N 2,X = A55(2)p2 

(A. 23) 

Pl= Ql.x 

= Mi,xx + ~2 Tl,x 

(A. 2 4) 


Matching the u displacement along the 1-2 interface implies 

Ui ( x -‘t)" U 2 ( x "T) 

hi h 

or Ui- T^bi = U 2 + * 2 ^b 2 


(A. 2 5) 


Combining the equations to eliminate the displacement and in- 


terlaminar stress terms leads to the following system of homogeneous 


coupled ordinary differential equations. 


Ni, x +N2,x » 0 

(A. 2 6) 

Mi.xx+aNi.^c-^M^O 

(A.27) 

M 2 ,x X - 2 2 N 2 ^-'^|M2=0 

(A. 28) 

Ni hiMi N 2 h 2 M 2 

Aii(i) ~ 2Dn(i) " Ai 1 ( 2 ) ' 2Dh(2) “ 

(A. 2 9) 


The solution is assumed of the form 


Ni] 

f Al 1 

n 2 

> _ J A 2 L 

Mi 

“1*3 

m 2 > 

IA 4 J 


r e sx 


(A. 30) 


Substitution of this solution into Eqs. A.26-A.29 leads to an eigenvalue 
problem with the following characteristic equation. 
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s[Bis 4 +B 2 s2+B 3 ] = 0 (A.31) 

where the B's involve the stiffness and thickness parameters A, D and 
h. For the material system and ply stacking sequence considered, 
B2 2 >4 BiB 3. Henee, the roots can be written as 


s=0,±' 


-B 2 ±VB2 2 -4BiB 3 

2Bi 


(A.32) 


Only the zero and positive roots of eq. A.32 are considered as they give 
solutions decaying exponentially from the crack tip. Then, the axial 


force and moment resultants can be written as 

Ni = aie s i x + a 2 eS 2 x + ai 
N 2 = -aie s i x - a 2 eS 2 x + a 2 
Mi = aikie s i x + a 2 k 2 e s 2 x 
M 2 = aik 3 e s i x + a 2 k 4 e s 2 x 


(A. 33) 
(A. 34) 
(A. 3 5) 
(A. 36) 


The k parameters in the above solutions involve the eigenvalues 
and the stiffness coefficients (A,D). For example, we have the defini- 
tion for ki as 


ki = 


hi o 

T Sl2 

^ 5 ( 1 ) 2 

Dn(i) _Sl 


(A. 3 7) 


Using the equilibrium Eqs. A. 10, A. 12 and A. 14 along with the 
applied axial force P and specimen width b, the axial force resultants 


and interlaminar stresses can be written as 

N 1 =aie s i* + a2e s 2x + ^ Aii ^l i|2) 

P Ah( 2 ) 

N 2 =-aieV - a 2 eS 2 * + jb A 11(1 )+An| 2 ) 
Ti = Ni, x = aisie s i x + a 2 s 2 e s 2 x 
Pi = Mi >xx +" 2 _ Ti (X 


(A. 38) 


(A. 39) 
(A. 40) 
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(A. 41) 


= ^ki+^jaiSi 2 e s i x +^k2+^ja2S2 2 e s 2 x 


hr 


The constitutive equations are used to write down the displace- 
ment solutions. The arbitrary constants associated with the displace- 
ments and rotations are determined from the matching conditions be- 
tween sublaminates 1 and 2 and the end conditions. Proceeding to 
sublaminate 3, the governing equations are 


N 3 ,x = 0 

(A. 42) 

Q3x + P 3 = 0 

(A. 43) 

s 

CO 

h 

I 

<D 

CO 

II 

O 

(A. 44) 

N 3 = Ai i(2)U3,x 

(A. 4 5) 

Q3 = A55(2)P3 

(A. 46) 

M 3 = Dh( 2 ) P 3 .X 

(A. 47) 


Matching U at the 2-3 interface and applying the boundary con- 
dition at the free end, N 3 (a) = 0 , gives 

N 3 = 0 

d ^ SiAn( 2 ) S 2 Ah( 2 ) d 

In order to solve for the bending moment, Eqs. A. 44, A. 46 
are combined to yield 

M 3 = cpi sinh ©3 x + (p 2 cosh ©3 x 
where ©3 is defined by 

_ a / ^55 [2 1 

0)3 " V D 1 1(2) 

Since the (3 matching conditon cannot be used at the 2-3 
the (remaining) boundary conditions are 

M 3 (a) = 0 
M 3 (0) = M 2 (0) 

The cps can be solved using the boundary conditions A. 52 and A. 53 as 


(A. 48) 
(A. 49) 

and A. 47 
(A. 50) 

(A.5 1) 
interface, 

(A. 52) 
(A. 53) 
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92 = aik3 + a2k4 (A. 54) 

91 = -92 coth (03a (A.55) 

The solution for sublaminate 3 can be completed by writing the ex- 
pressions for Q3, P3 and P3 based on the M3 solution. 


The equilibrium equations for sublaminate 4 are 

N 4 , x = 0 (A. 56) 

Q 4t x = 0 (A. 57) 

M 4tX - Q 4 = 0 (A. 58) 

The constitutive relations take the form 

N 4 =An(i)U 4 ,x (A. 59) 

Q 4 = A55(i)(p 4 +W 4tX ) (A. 60) 

M 4 = Dh(i)0 4i x (A.61) 

P 

Using Eq. A.56 with the boundary condition N 4 (a)=2^ yields 

N 4 = (A. 62) 


Similarly, using Eq. A. 57 with Q 4 (a)=0 results in 

Q 4 = 0 (A. 63) 

Matching Mi and M 4 at the 1-4 interface and using Eq. A. 58 


gives 

M 4 = aiki + a2k2 (A.64) 

The U 4 displacement is obtained by integrating Eq. A. 59 and 
using the displacement matching boundary condition U 4 (0) = Ui(0). 

°*=Ans(^ x+ fi + if) +a 3 (A - 65) 

Similarly, integrating Eq. A.64 and setting 0 4 (a) to zero gives the 
solution for 0 4 . Using the solutions for Q 4 and 0 4 and the boundary 
condition W 4 (0)=0 in Eq. A.63 yields the solution for W 4 . 
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In order to determine ai, a2 and a3, the following boundary 


conditions are used. 

Ni(0) = ^ 

Pl(0) = p 4 (0) 
Ui(-L+a) = 0 


It is convenient to define the following parameters. 



(A. 66) 

02 = Dnoi 

(A. 6 7) 

03 = A55UI ( k2 + 2 ) 

(A.68) 

9 « = Dn (1) 

(A. 69) 

0d = 03 • 01 + (04 - 02)a 

(A. 70) 

The nominal (far field) strain is given by 


P 1 

£ “ 2b An(i)+ Ah(2) 

(A. 71) 

The a parameters are obtained as 


. 03+04a 

a l =Aii(2)£ 

(A. 72) 

, 0i+02a 

a2 = -An(2)£ 

(A. 73) 


as = 6 (L ' a) • sIaTTo) e ' Sl(L_a) 


S2A1KD 


e -s 2 (L-a) 


(A. 74) 


The specimen compliance C is defined as the ratio of specimen 
extension to applied load. This is obtained as 


2U 4 (a) 
= p 
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a l 3.2 
+ 77 + ^ + a 3 A na) 


(A. 7 5) 


~ PAn(i) [2b si S 2 
The total energy release rate associated with the crack 
(delamination) growth under a constant load P is given by 

(A. 76) 


GT_ 2b da 


Substituting the compliance expression from Eq. A. 7 5 in Eq. 
A. 76 yields the following expression for the total energy release rate. 

P 2 f 1 1 



Ah(i)+Aii(2) 


+Ii 


•bl 


where 


0203-0104 rl-e-si^ 3 ) l-e-^tL-a) \ 

U= X 0 d 2 J 


I2=X 


with 


(03+04a)e~ s it L ~ a ) - (0i+02a)e~ s 2^ L ~ a ) 
0d 


Y = _ __I . A ll(2) 
K Aii(i)+Aii(2) And) 


(A. 77) 

(A. 78) 
(A. 79) 

(A. 80) 


The individual fracture mode contributions to the energy release 
rate can be calculated using the virtual crack closure method, based on 
the interlaminar stresses and displacements in the vicinity of the 
crack tip. From the assumed plane strain condition, the mode III con- 
tribution is zero (Gm=0). The mode II energy release rate, Gn, is cal- 
culated using the virtual crack closure technique while Gi is evaluated 
using 

Gi = Gt - G 11 (A. 81) 

Gn is calculated from the interlaminar shear stress and relative 
sliding displacement as 


lim 1 f 
G n = 8->0 jTifr-SJAnMdx 


(A. 82) 
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* 


In the absence of a singularity in the stress field, the limiting 
process leads to the trivial result Gn= 0 . Hence, the limit is calculated 
as 8 tends to some finite value, say A. The value of A is chosen de- 
pending on the decay length associated with the problem i.e. the 
length within which the presence of the crack significantly alters the 
specimen response in comparison with the corresponding far field 
values. The decay length in this problem is dependent on the eigen- 
values si and S2- The following value of A has been chosen in order to 
reasonably fulfil the decay length criterion. 



(A. 83 ) 


The relative sliding displacement Au is based only on the differ- 


ence U4-U3 so that the kinematic condition of zero relative displace- 


ment at the crack tip is fulfilled. This also simplifies the calculations. 
If the true value of Au (based on U4-U3) is used, the (3 mismatch at the 
3-4 interface leads to a kinematically inadmissible displacement dis- 
continuity at the interface. This discontinuity causes a non trivial lim- 
iting value G11 as 8-»0. But this value is an artifact of the modeling as- 
sumptions and cannot be used as the true value of Gn. The mode II 
energy release rate component, using Au=U4-U3, is obtained as 

Gn = ^ (A. 84 ) 

where the parameter I3 depends on Aim), An(2), si, S2, ai, a2, A and 
the specimen nominal strain e. 
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Table I Summary of Results 


number of 
90° plies 

Gt 

J/m 2 

Gu/Gj 

<?c 

MPa 

e c 

% 

1/2 

2.404 

0.276 

1313.9 

1.6747 

1 

6.752 

0.275 

784.0 

1.1685 

2 

22.849 

0.267 

426.2 

0.8058 

3 

51.049 

0.261 

285.1 

0.6427 

4 

93.603 

0.256 

210.6 

0.5444 

6 

228.871 

0.250 

134.7 

0.4264 

8 

440.065 

0.247 

97.1 

0.3555 



Fig. 1 Specimen Cross Section 
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Fig. 2 Modelled Regior 



Fig. 3 Generic 
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Fig. 4 Sublaminate Forces and Coordinate Systems 








Fig. 7 Interlaminar Shear Stress Distribution 



Fig. 8 Interlaminar Normal Stress (Peel Stress) Distribution 
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Fig. 9 Energy Release Rate Comparison 



Fig. 10 Critical Delamination Growth Strain Variation 




